arXiv:1508.01276vl [cond-mat.stat-mech] 6 Aug 2015 
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A rigorous derivation of nonequilibrium entropy production via the path-integral formalism is 
presented. Entropy production is defined as the entropy change piled in a heat reservoir as a result 
of a nonequilibrium thermodynamic process. It is a central quantity by which various forms of 
the fluctuation theorem are obtained. The two kinds of the stochastic dynamics are investigated: 
the Langevin dynamics for an even-parity state and the Brownian motion of a single particle. 
Mathematical ambiguities in deriving the functional form of the entropy production, which depends 
on path in state space, are clarified by using a rigorous quantum mechanical approach. 
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INTRODUCTION 

The fluctuation theorem (FT) is an important princi¬ 
ple in nonequilibrium statistical mechanics. It was first 
discovered two decades ago for a deterministic molecular 
dynamics [1], but was later proven to hold for stochas¬ 
tic dynamics ns. It deals with thermal fluctuations of 
time-accumulated quantities such as work, heat, and en¬ 
tropy production (EP). In the usual statistical-mechanics 
approach, only average quantities are of interest because 
fluctuations are negligible in the thermodynamic limit. 
However, owing to modern technologies, one can inves¬ 
tigate the thermal motion of a small system on a nano¬ 
scale that shows large fluctuations. For example, one can 
observe that the thermodynamic law is violated in rare 
events occurring due to extreme fluctuations. The FT 
provides a rigorous relation for the probability distribu¬ 
tion function of such a fluctuating quantity and has been 
proven for various small-system experiments 

The path-integral theory for the thermal motion was 
developed by Onsager and Machlup in 1953 [16J and has 
been used extensively to study the fluctuation theorem 
and related issues. Finding the EP from the path prob¬ 
ability is essential in order to prove the FT. Consider a 
particle evolving along a path q(f) in state space under a 
time-dependent protocol A (t) for time period 0 < t < r, 
starting from q 0 and going to q T . The state vector q 
may be the position x in an overdamped case or the 
position-momentum pair (x, p) in an underdamped case. 
Infinitely many paths connect the initial and the final 
states, among which a single path is determined stochas¬ 
tically from the given dynamics. Given a path q(f), the 
time-reverse path is given by 


q 0) = 


be that for the time-reverse path under the time-reverse 
protocol A(f) = A(r — t). Then, the EP is given by 


Ai'env — k B In 


H[q(£); A(f)] 
H[q(t); A(t)] 


( 2 ) 


where ASenv is another name for the EP denoting the 
entropy change in the environment (heat reservoir), and 
ks is the Boltzmann constant. This relation is called the 
detailed fluctuation relation [H in the literature, which 
was first derived by Schnackenberg for the master equa¬ 
tion for discrete states 


17 ]- 


Then, the change in the total entropy of the system 
and the heat reservoir is given by 


A5 tot = k B In 


/o(qo)n[q(t); A(t)] 
p(q-r)n[q(t); A(f)] 


( 3 ) 


where p(qo) are p(q r ) are the initial and the final proba¬ 
bility distribution functions (PDF), respectively. Rec¬ 
ognizing the system entropy as the Shannon entropy, 
— In p(x) 0 , the change in the system entropy is given 
by AS'sys = fc B [lnp(qo) - p(q T )], so A S tot = AS sys + 
ASenv ■ Then, the FT for the total entropy change can 
be easily proven to be 


(e fc s lAS tot) 


dqop(qo) J £>[q(i)]n[q(f); \{t)] 
p(q T )n[q(t); A(t)] 


( 4 ) 


x(t — t ) (even parity) 

(x(r — t), — p(r — t)) (mixture of even and odd 

( 1 ) 

where x has even parity and the momentum p has odd 
parity in time reversal t —> r — t. Let n[q(f); A(f)] be the 
conditional probability for the system to evolve along the 
forward (given) path under the protocol and n[q(f); A(f)] 


p(qo)n[q(t); A(f)] 

= J dq 0 p(q T ) Jornwimm = 1 , 

where the Jacobian for the change of variable q —> q 
is equal to 1 and p(q r ) is the initial PDF for the re¬ 
verse path. Using Jensen’ s inequality, one can get 
ot) — 0) which confirms the second law of ther¬ 
modynamics. Other types of FT’s can also be proven 
in similar ways where the use of the detailed fluctuation 
relation in Eq. 0 is crucial [19]. 

In this paper, the detailed fluctuation relations for two 
kinds of the stochastic dynamics are derived: (i) the 
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Langevin dynamics for an even-parity state and (ii) the 
Brownian motion of a single particle. The dynamics (i) 
is described by 

x = g(x, X(t),t) + £(t) , (5) 

where £ is the noise vector with zero mean and with 
the variance given by (£(t)£(tf)) = 2D S(t — t'). The 
outer product cd for vectors c and d denotes a matrix 
with element c m d n for the rrfi 1 row and the rfi 1 column. 
D is called the diffusion matrix, and is symmetric and 
positive-definite, g is the drift term, is proportional to 
the applied force in the overdamped limit, and depends 
on the protocol as well as the state of the system. The 
dynamics (ii) is described by 


x = p/m 
p = —G • p/m - 


f(x,p, X(t),t) + C(t) 


( 6 ) 


where x and p are the position and momentum of the 
particle, G the friction matrix, and m the mass of the 
particle, f is the applied force on the system, which 
generally depends on the position and the momentum 
as well as the protocol. S is the diffusion matrix for the 
noise vector £, given from (£(t)£(t')) = 2S 8(t — t'). In the 
overdamped limit with large G or equivalently in the zero- 
mass limit, (ii) is found to reduce to (i) if f is independent 
of p. In reality, the momentum-dependence of f can be 
found in active matters [1301 and the feedback con¬ 
trol [23-1251. Many unprecedented features have recently 


been found for such an odd-parity force [iH0 . but not 
for any even-parity force. Only non-multiplicative noise 
where the diffusion matrix is independent of x or p is 
considered, and a study of the multiplicative case will be 
presented in a separate place. 


LANGEVIN DYNAMICS FOR AN 
EVEN-PARITY STATE 

The differential equation for the PDF associated with 
the Langevin equation, Eq. ©, is given by the Fokker- 
Planck equation 0 |: 

dtP(x,t) = -H F pp(x,t) , (7) 

where <9 t = d/dt. Hp p is a non-Hermitian operator given 

by 

Hp P = dx ■ [g - D • <9 X ] , (8) 

where <9 X = V x . Introducing an operator a = * _1 9 x for 
i = y/—l conjugate to x, one can write 

Hp p = ih • g + a • D • a . (9) 

Using the quantum-mechanics formalism, p(x, t ) = 

(x| p{t)) = (x|e tH FP |x 0 )(xo|p(0)), where the integra¬ 


tion over Xo is implicit following the Einstein conven¬ 
tion. The propagator (x|e tH FP |xo) in the quantum- 
mechanical sense turns out be the conditional probabil¬ 
ity P(x,t|xo,0) where stochasticity (normalizability) is 
guaranteed by / dx(x|iLpp |xo) = 0, which can be seen 
easily from the total derivative in Eq. ©. 

For 0 < t < r, let us discretize time by by using the 
interval A t = t/N as tj = jAt for integer j. For 0 < j < 
N , to = 0, fjv = t. In the TV — > 00 (At — > 0) limit, one 
can write 


P(x t ,t|x 0 ,0) = 


N—l 

n <ix i 

3= 1 


N 

n 

3= 1 


(x i |e- AtH FP| Xj _ 1 ) , 


( 10 ) 

where x ; = x(tj). A set {xo, xi,..., Xj,..., xjv} is a dis¬ 
crete representation of a path x(t). The integration over 
all intermediate states x ? for 1 < j < N — 1 corresponds 
to the integration over all paths connecting xo and x r . 
Then, one can identify the path probability, II[x(t); A (t)], 
as nf = i(x,|e A<H FP|xj_i), where A (t) takes the value 
Xj_i at time T ? _i for the interval between tj_i and tj. 

One can change the order of a in TFpp as Ta • g = 
ig ■ a + 9 X • g, by using the commutator [a, g] = * _1 9 x g. 
Changing the order partially, one can write 


-®FP = (1 — a )Hpp + aHpp , (11) 


where 0 < a < 1 and Ppp = ig ■ a + d x ■ g + a • D • a. 
For eigenstate |a) of a, 1 = |a)(a| and (x|a) = (a|x)* = 
e ! a - x /( 2 7r )d /2 j n d-dimensional space. Then, one can find 
the expression in detail for the propagator between tj -1 
and tj as 


(xj|e AtH FP |xj_i) 

- ( x j l a )( a l [1 - Ai (! - oi)H Fp ] |xj_i) 

— ( x j | aAtHpp | a) (a| Xj _ 1 ) 

~ f ^ a ^ — A^a-D-a-f-za-(x,— iAta-g^ — Ata(d x -g) i 

~ J ( 2 n) d 


[det(47rAtD)] 1 / 2 


e-^h^bD-Lh^-AtaOx-g) 


( 12 ) 


where = (xj — Xj-\)/At — g* a) . g^ is an interme¬ 
diate value between g at Xj_i and x, and is defined by 
(l-a)g i _i+ag J for g k = g(x fc , Aj_i) for both k= j -1 
and j. Special cases are a = 0 (prepoint), a = 1/2 
(midpoint), and a = 1 (postpoint). The meanings of the 
terms in the brackets are clear. Note that an arbitrary 
use of a for the path integral leads to the same result up 
to higher-order correction O(At). 

Then, the path probability is found as 


n [x(f); A(t)] 

= lim Me-^= lAt [^ ) D " h ^ +a{d ^ 

N—too 

— _^ e -/ 0 T dt[ih(x,x,i).D -1 -h(x.x,t)+a9 x .g] 


( 13 ) 
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where J\f = [det( 47 rAtD)] _JV / 2 and h(x, x, t) = x— g. The 
last line is written in the continuous-time limit while 
it is only well-defined in the discrete-time representa¬ 
tion. The integration of Eq. m over all paths and over 
the initial state, given the initial PDF, goes to 1, which 
can be written as 1 = f dx 0 p(x 0 ) J D[x(f)]II[x(f); A(f)]. 
This normalization property is used in Eq. ©, where 
_D[x(f)] = limjv->oo-A/ - IljLi dx.j • For the prepoint rep¬ 
resentation, the continuum limit can be seen to be the 
same as the noise distribution by noting x — g = £(i), 
agreeing with common intuition. The seeming ambiguity 
of the path integral, particularly in treating the diver¬ 
gence term d x ■ g, can be well resolved by a correct use of 
the a-representation in Eq. urn the choice of which is a 
matter of taste. The path integral can be exactly carried 
out when g is linear in x, which was done for the linear 
diffusion problem using the prepoint representation [ijoj ■ 
It was also done using the midpoint representation for a 
similar linear problem of the injected and the dissipated 
powers 


31|. 


The time-reverse path is given by {0 < j < JVjxj = 
xjv-j}. For the time interval between tj-\ and fj, the 
drift term is denoted as g*, = g(x£,Aj) for both k = 
j — 1 and j, where the protocol has a fixed value at the 
later time, opposite to the forward path. This is the 
strategy for the dynamics to have a common protocol 
for an overlapping segment of path. Replacing j with 
N — j + 1, the transition occurs reversely from Xj to Xj_i 
under a fixed protocol value Aj_i. The corresponding 
forward transition occurs from Xj_i to Xj under the same 
protocol. Then, the probability for the reverse path can 
be written by using a different discretization-parameter 
7 . One find that 


n [x(t);A(t)] 


N—> 00 


= lim A 7) ' d 1 ' h f 7 ) +7(&c-gb-i] 


N—> 00 


(14) 


= Y Ax J • D 1 ' gf /2) > (!5) 

i =1 

where Axj = Xj — Xj_i, and higher-order corrections of 
order (At ) 1 / 2 are neglected. Note that g(“^ ~ gi 1 / 2 ) 4 - 
(a — l/2)(Axj • <9 x )g j- Then, one can show Axj • D ” 1 • 
(Axj • d x )gj = 2Tr[(9 x gj)D _ 1 (AxjAxj)], which goes to 
2A td x ■ g j when using AxjAxj = 2AfD. Therefore, the 
term with the factor (a — 7) vanishes. Note that the 
EP is uniquely defined by the midpoint representation, 
irrespective of the representations used for the forward- 
and the reverse-path probabilities. The importance of 
the midpoint representation for the EP lies in the fact 
that gj 1 ^ is the Taylor expansion of g j up to O(Axj), 

so A Xj ■ D -1 • gj 1 ^ 2 "* can be found up to the second order 
in A Xj, which is the desired order O(At). For a simple 
example in one dimension, consider g = —kxi. Then, 
AiSenv is proportional to (Xj — Xj-i)kx^ 2 ^ = k(xj — 
Xj-i)(xj + Xj- 1)/2 = A;(a ; 2 — x 2 _ 1 )/2, which becomes 
0{At). 

In the continuous-time limit, one can write dSenv/dt = 
ks x ■ D _1 • g. The Fokker-Planck equation, Eq. m , can 
be written as dtp = — d x • j, with the probability current 
j = (g— D -d x )p. Let the steady-state PDF be pss oc e ~^ 
with a potential 4>(x). Then, the steady-state current is 
given by jss = (g + /3D • 9 x ^)pss- The detailed balance 
(DB) was found to be characteristic of equilibrium and 
to be measured by jss 0, 52, 53] • The DB holds if g = 
—/3D- d x (j) while it breaks if D ^ 1 • g has a non-conservative 
part fnc- Generally, in the presence of a time-dependent 
protocol, the drift term can be written as 

g = • [—d x cj>(x, A(t)) + fnc( x , t)] , (16) 

where the time-dependent protocol is supposed to change 
the potential; otherwise it can be included in gnc- 
(/fD ) -1 -g can be said to play the role of the force exerted 
on a particle in motion. From this point of view, one can 
write the generalized Clausius law as 


where h) 7 '* is given by h]' 11 in Eq. (fl2l) with the replace¬ 
ment of Xj by Xj, and hj 1 7 ' ) = — (xj— Xj_{)/At— gj 1 1 \ 
Using Eqs. m and m, one can find the expression 
for the EP in Eq. © as 


AS'env 


N 


= £ 


3 =1 

N 


A xj • D" 1 • 


1 


g ■ Q) +g - 1 7) 


- (a - 7 )A td x • gj 


Y, [Ax.7 • D 1 -g 


-1 . „(V2) 

'3 


3 =1 


~{a - 7 ) ( A td x ■ gj - - A Xj • D 1 • (A Xj ■ <9 x )gj 


dSe nv 1 dq> • 
dt ~ T ~~dt + 


Q 

T ’ 


(17) 


where d(f>/dt = d x <f> ■ x + {d(j)/d A)A has been used. W is 
interpreted as the rate of work produced by the nonequi¬ 
librium sources, and is given by 

W = fnc • x + —^ A . (18) 

a A 

Q is interpreted as the rate of heat produced in the 
heat reservoir, which implies the generalized first law of 
thermodynamics, d(j)/dt = W — Q. One can envisage a 
stochastic process governed by the Langevin equation as 
a thermodynamic process, which maybe is out of equilib¬ 
rium. 
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In the overdamped limit for the Brownian motion of 
a particle described in (ii) in the Introduction, one can 
neglect the inertia part p in Eq. ([ 6 ]) and finds x = G -1 -f+ 
G _1 • Then, one finds g = G -1 • f and D = G _ 1 S(G _1 ) t . 
For the heat reservoir in equilibrium at temperature T, 
the Einstein relation S = /3 -1 G holds so that one can find 
x • D” 1 • g = /?x • f . Identifying f as — 9 x </>(x, \(t)) + fnc, 
the above result in Eq. 03 has a physical meaning. 


BROWNIAN MOTION OF A SINGLE PARTICLE 

The differential equation for the PDF associated with 
Eq. ([ 6 ]) is given by the Kramers equation [29|, which is 
the Fokker Planck equation for p(x,p): 

<9 t/ o(x, p, t) = -H K p(x, p. t) (19) 

with 

tf K = d X -- + V (-G--+f-S-d p ) , (20) 

where d p = V p . Using the quantum-mechanical formal¬ 
ism, it can be rewritten as 

H k = ik ■ — + ib • f-G • 3. + f) + b • S • b , (21) 

where a = i~ 1 d x is the same operator as in the last 

section and b = i~ 1 d p is the conjugate operator for p. 
Using the commutation relation, [b, f] = i _ 1 <9 p f, one can 
write 

H K = (1 - a)H K +aH K , (22) 

where 

H k = i— -a + if-— ■G t +f > ) -b + b-S-b 

m \ m J 

+dp-f. (23) 

to 

Let |ab) be a simultaneous eigenstate of a and b that has 
the eigenfunction given by (xp|ab) = (27r) _d e ia x+ * b ' p . 

In the same time-discretization scheme as in the last 
section, the path probability can be found for a discrete 
path {0 < j < iVjxj,pj}. Using the a-representation, 
one can find 

N 

n [x(£), p(t); X(t)\ = lim [|(x J p i |e _Afff K|x j _ip :; _i) 

N—too 

3 =1 
N 

= Jim n [( x jPjl ab )( ab l [! - Ai (i - a ) H K\ l x I-i x i-i) 

—>oo A A L 
3 =1 

- ( x j Pj I °A/// k | ab) (ab |xj_ lPj -_i) 
jv . , 

— lim TT / __f_ P ia [( x i- x i-i)- Ai p)“ ) /"i] 

- (2n) d 


^ b ^ — Aib-S-b+zb-(pj— pj_i )—— Gp(“^/m+f) a ^) 


xe 


(27r) d 
Ata[(TrG)/ m—d p -fj] 


lim TT 5 ^ j ~ Xj_1 ~~ 

AF —-A- 


[det( 47 rAtS )] 1 / 2 


-AtEf=i[r'“ ) -S- 1 -r(.“V4-a(TrG/m-a p .fU] 


N- 

3 =1 

xe ““^ 3 =i 1 ‘3 


(24) 


where = (pj — Pj_i)/At +G • p^/m — fj a \ In this 

expression, fj 0 -* = (1 — a)fj-i + afj is an intermediate 
value of force f*, = f(x*,,p*,, Aj_i) for k = j — 1 and j, 
where the protocol has the fixed value at the earlier time. 
In the continuous-time limit, Eq. AMD can be written as 

II ( <5(X - p/ m )e- fo dtir-S-^r/A-aCTl-G/m-d^f)] up tQ the 

normalization constant. 

The corresponding time-reverse path is given by {0 < 
j < IV|Xj = xjv-j, Pj = —p N—j }■ For a time inter¬ 
val between tj-i and tj , the force is denoted by f*. = 
f(xfe, pfc, A j) for k = j — 1 and j, where the protocol has 
the later time. Replacing j by N — j + 1, the transition 
occurs reversely from (xj,—pj) to (xj_i,— Pj-i) under 
a fixed protocol Xj-i- The corresponding forward transi¬ 
tion occurs from (xj_i,pj_i) to (xj, pj) under the same 
protocol. The probability for the time-reverse path can 
be found in 7 -representation as 


n[x(t),p(£); \{t)] 

_ tt H^N-j+1 — xjv— j — Atp N _- + 1 /m) 

at- 3 -oo ii [det( 47 rAtS )] 1 / 2 

xe~ At S^= 1 [ ? w-j+ 1 ' s_1 ' f w-i+i/ 4_ T'( TrG / m_a p'^-i+i)] 

u A S (*j - x j-i - Atpf^/m) 
[det( 47 rAtS )] 1 / 2 

xe -AiEf=i^ ) -s- 1 -?^V4-7(TrG/m-a P -L'- 1 )] , ( 25 \ 

where r ^ is given by r^ 7 -* with the replacement: 
Xj,pj —>Xj,pj. Then one finds = (P iJ —Pj-i)/A£ — 
G-p j 1-7) / to - for fj = f(Xj, —pj,Xj). Here, 

p^l\ = p * 1 7> has been used. 

It is useful to write f iev = (f+f)/2 and f br = (f—f)/2. 
Under p —> — p, the former is symmetric while the latter 
antisymmetric. Write Ap j = p j — Pj-i- Then, the EP 
in Eq. m can be written as 


fcg 1 AS'env 


N—too 

E 

3 = 1 


(Ap 


, - Aff; 


xev(a) 


•S' 


+ (A P ,-Atf; ev( 1 - 7 ) ).s- 1 . \G. 
-Ap, ■ S” 1 ■ (ff v(a) - 





TO 


» 


- f 


irr(a) 


TO 


- f; 


irr(i-7) 











5 


-2 At 


a — 7 


TrG — d p ■ (afj + qfj) 


One should keep the terms in the bracket up to O(At) 
and note that the outer product ApjAp,, is equivalent 
to 2A/S. One can expand fj“^ = f j 1 T) + {a — l/ 2 )(Ap,- • 
<9 p )fj and do the same for the other vectors. After rear¬ 
ranging the terms, one can simplify the EP in the time 
interval between tj- i and tj as 

( 1 / 2 ) 

k^AS pnv i = -Ap, • S _1 G • ^-+ Atf, • S _1 G • Si 

+Ap, • S " 1 • f irr(1/2) - Atf irr • S _1 G • Si 

J J m 

— AifJ ev • S - 1 • f] rr - A td p • ff v 
+remainder , (27) 

where the term ‘remainder’ is given by 

(a - 7 ) (a t^- - ^A P j • S _1 G • (Ap, • <9 P )^ 

+ ^ A Pj • S_1 • ( A Pj • d P ) ((« - 7)ff r + (a + 7) f J eV ) 
-AtOp- ((a + 7)fi' ev + (a-7)f] rr ) , 


which can be shown to vanish by using Apj • S _1 G • (Ap, ■ 
9 p )p,- = 2TrG/m and Apj • S _1 • (Ap, • <9 p )f, = 29 p ■ fj. 
Therefore, the EP is uniquely defined from the midpoint 
representation, independent of the specific representa¬ 
tions used for the path probabilities. Note that only the 
terms coupled with Ap„- should be written in the mid¬ 
point representation. 

In the continuous-time limit, the EP can be written as 


doe nv 

dt 


TrQT - 1 


d5an 

dt 


(28) 


The matrix T -1 of the reciprocal temperature is defined 
as /cbS - 1 G. The matrix Q of the heat production rate is 
defined as — (p — f)p/m, the trace of which is the real 
rate of heat produced in the heat reservoir. The second 
term is the rate of the anomalous EP, which was found 
in a recent work [28|. One can get 


TrQT - 1 

dSan 

dt 


— (p ~ f) ■ S -1 G • — , (29) 

m 

f irr . s -1 . (p _ G . -B. — f rev ) - dp • f r fyo.) 


fixed temperature T, has been of much interest [ISEDI. 
In this case, one has G = yl and T = T I, and finds 
S = 7 fcsTI, which is known as the Einstein relation or 
the fluctuation-dissipation relation. Therefore, one can 
find TrQT -1 = — (p — f) • p/mT = Q/T , where Q is the 
rate of work extracted by the force, — yp/m + £, exerted 
by the reservoir, which is, by definition, the heat produc¬ 
tion rate. One can write f = — <9 x E(x, A (t)) + fnc( x , p)- 
In the discrete-time representation, one can find 


-(p-f)-- 

m 

1 A Pi 


A 1/2) + (_a xl y. + f „ .). P L 

m 


m At 3 

P/ P/~i dV dv { Pj 

2mAt dt dX ncj m 

-§+*-<>■ 


(31) 


where E = p 2 /(2m) + V (x, A) is the total energy and W 
is the work production rate, which is composed of two 
contributions: 


(J V • 

w = —-A + f nc -X . (32) 

d A 

The first term was recognized by Jarzynski Q, and the 
second term is the rate of work done by the nonconser¬ 
vative force. The final line in Eq. m is consistent with 
the first law of thermodynamics. If there is no time- 
dependent protocol and no nonconservative force, only 
heat is dissipated by the amount of energy change. In 
the steady state, therefore, no heat is produced, which is 
characteristics of equilibrium. On the other hand, if work 
production occurs from a time-dependent protocol or a 
nonconservative force, persistent productions of work and 
heat can occur even in the steady state with (dE/dt) = 0, 
which is a fundamental difference between an equilib¬ 
rium steady state and a nonequilibrium steady state. In 
the absence of a momentum-dependent force, the Clau¬ 
sius law is, indeed, satisfied, dSeuv/dt = Q/T. How¬ 
ever, in the presence of a momentum-dependent force, an 
anomalous production ASan added to Q/T appears. An 
anomalous feature is even more serious if an irreversible 
part f 111 exists, which can be seen in Eq. m- There is 
yet no clear understanding of ASan beyond its existence. 


CONCLUSIONS 


Here, the midpoint representation is implicit while it 
should be explicit for a rigorous calculation of the EP. 

Equation (l29l) gives the expression of the EP for a 
general heat reservoir that might not be in equilibrium 
for asymmetric T _1 = ks S^ 1 G. The Brownian mo¬ 
tion of a colloidal particle in a liquid, where the liq¬ 
uid plays the role of an equilibrium heat reservoir at a 


A rigorous expression for the EP, the detailed fluctu¬ 
ation relation, that is composed of important thermo¬ 
dynamic quantities such as energy, work, and heat, is 
derived. The result is shown to be uniquely given by the 
midpoint representation independent of any representa¬ 
tions used for the forward- and the reverse-path proba¬ 
bilities. 
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The exact rule for the discrete-time representation is 
very important in the calculation via the path-integral. 
In order to find the average value of a thermodynamic 
quantity, one needs to carry out the path integral with 
the path probability. Then, one can use any representa¬ 
tion for the path probability provided in Eq. m or m- 
If the quantity to be averaged is path-dependent, such as 
EP, work, and heat, it should be expressed in the mid¬ 
point representation. For example, Ap,, • p"' 1 /m goes 
to the correct kinetic energy change (p^ — p|_ 1 )/(2m) 
only if a = 1/2 is used. If a different value of a is 
used, the correction is (a — l/2)Ap j ■ (Apj • <9 p )pj/m = 
2(a — l/2)Af(TrS)/m, which is not negligible. The sec¬ 
ond law of thermodynamics can be expressed by using 
the positivity of (S'tot) =1 dxdpp(x, p) -1 j lrr • S _1 • j lrr , 
where j^ 11 = (— G-p/m+f^ 11 — S-d p )p. The midpoint rep¬ 
resentation of the EP is necessary for the derivation |28| . 

The detailed fluctuation relation in Eq. © is an in¬ 
gredient to the FT for the total entropy production 
shown in Eq. ©. One can also use it to prove other 
FT’s, for example, for work production and heat pro¬ 
duction. There are two types of the FT. One is called 
the integral FT (IFT) of type ( e~ R ) = 1 for a time- 
accumulated quantity R. The other is called the detailed 
FT (DFT) of type P F (R)/P R {-R) = e' iR . The former 
is a necessary, but not a sufficient, condition for the lat¬ 
ter. The IFT can be proven in a way similar to that 
used for Eq. ©• The DFT can also be proven by us¬ 
ing -Pp(-R) = (S(R — R[q(t), A(t)])}, where the bracket 
denotes the integral over all paths and the initial state. 
P r (—R) is given similarly by (S(R+R[q(t), X(t)})), where 
the subscript R stands for the process under time-reverse 
protocol. See Ref. [34| for the detail. 

The investigation in this study is restricted to stochas¬ 
tic motion due to additive noise. A derivation of the 
detailed fluctuation relation for multiplicative noise will 
be presented in the near future. 

This work was supported by the research fund of My- 
ongji University, 2013. 
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